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Antenna array of large scale have been examined for different applications 
including 5G technology. To get better data rate or a reliable link substantial 
number of antenna arrays have been utilized to provide high multiplexing 
gains as well as array gains with high directivity. In this paper a simple but 
efficient implementation technique of using sub-arrays for the improvement 
of large-sized uniform arrays. By repeating a small sub-array multiple times 
large arrays can be designed. This implication of utilizing small array 
simplifies the design of a larger array which allows the designer to 
concentrate on the smaller sub-array before assembling larger arrays. So, by 
investigating the sub arrays the performance and radiation characteristics of 
large arrays can be anticipated. The array-factor for a planar sub-array of 2x2 
(4 elements) is analyzed using Mat-lab software and then a large array is 
formed by placing the 2x2 sub-array in different configurations in a 
rectangular arrangements up to 8x8 planar array. And then the results are 
validated with CST (Computer simulation technology) simulation results. In 
this way the array-factors, directivities, HPBWs, and side lobes of the 
constructed large arrays are analyzed and associated with the small sub-array. 
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1. INTRODUCTION 

Large-Scale Antenna Array (also known as Massive MIMO or Hyper MIMO) system is a 
vital notion in wireless communications research areas, which highlights the importance of advanced 
5G technological network architectures. At the base stations a huge number of array antennas operate 
consistently in order to direct the signals into smaller regions of space. Additionally, two or more antennas 
is outfitted with the consumer gadgets to facilitate the reception of signals, consequently there is 
a considerable enhancement in electricity effectivity as well as the throughput [1-4]. As large antenna arrays 
generally resonate in millimeter wave frequency ranges, their position can be confined in a very small area 
with a variety of array configurations i.e. planar, circular, spherical or cylindrical [5-6]. 

The configuration of the array placement, the range of elements, the characteristics of the radiation 
pattern, the inter component spacing and mutual coupling amongst the array factors each make contributions 
in the ultimate performance of any massive MIMO system. The consistent frequencies that are used for 
Massive MIMO structures are 3-6GHz, 27-28GHz and 60-70GHz, whilst the kind of antennas investigated 
can be shaped as patch antennas, horn antennas and also dipoles in some cases. In Figure 1, a 2x2 rectangular 
antenna array configuration is used in cell phones as shown. Large planar arrays such as 4x4, 8x8, 16x16, and 
32x32 can additionally be geared upinto small areas of less than one square centimeter at the base stations. 
Planar arrays are characterized with high gains and very directive beams in the preferred directions whilst 


Journal homepage: http://beei.org/index.php/EEI 








Bulletin of Electr Eng and Inf 


ISSN: 2302-9285 


□ 907 


displaying low side lobes in undesired directions. Moreover, as the gain and directivity of an antenna are 
related to each-other, the array factor can be effortlessly calculated from that relation [7-9]. 



Figure 1. Array synthesis when using sub-array concept [9] 


There are quite a few techniques amongst the researchers have been improvised in order to study the 
characteristics and behavior of massive antenna arrays radiation patterns. One of the most famous strategies 
of them is used in this venture is the notion of sub-arrays. Sub-arrays are regarded as factors in a large array. 
This provides an easier analysis as the designer will focus on the small sub-array design and construct large 
arrays by repeating those sub-arrays throughout the entire large array. According to [10-12], the overall 
pattern of an array arises from the multiplication of individual element pattern and the array factor. When 
sub-arrays are involved, the overall pattern can be achieved from the product of element pattern, the sub¬ 
array factor, and the full array factor determined by the amplitude of excitation, spacing between the 
elements as well as phasing between sub-arrays. 

Sub-array implementation is a very beneficial technique for designing massive array antennas. 
It simplifies the design by allowing the large array designer to concentrate on the smaller subarray design 
first and finally use the sub-array element as a part of any other larger array. Without indulging into full array 
study, this new approach implements the sub-array to construct a larger array. This process reduces large 
array problems into small array solution which has the potential to save a lot of CPU time and requirement of 
heavy memory as well. Generally, radiation characteristics of any large array can be easily predicted if the 
mutual coupling is ignored between the elements of array antenna. By implementing pattern multiplication 
approach of the array factor as well as the radiation pattern of single element would be adequate to forecast 
the array radiation performance and characteristics [13-15]. 

There are many other widely used numerical approaches available for large array analysis. Some of 
the analysis are based on infinite-array approach as well as the periodicity of the structure while others 
depend on the actual finite size array. Infinite array method ignores the edge effects while the simulation of 
the second approach is highly computational, complex and expensive when it comes to large array size. 
Research mentioned in [9] validates that subarrays can be used as a building block to facilitate the design of 
large antenna arrays. By controlling the subarray, the overall radiation pattern can be easily predicted. 
This can be achieved by multiplying the element pattern, sub-array factor and the overall array factor. Figure 
1 represents the levels of array factor synthesis. 

Any changes that may take place regarding the pattern of the elements of the array factor, due to that 
change subarray or full array factor will also be hugely affected in the final outcome of the array pattern. 
Changes can appear in the form of maximum (grating) lobes, variation of side-lobe levels or any other 
change into symmetrical dysfunction of the array pattern. Grating lobes appear due to distortions in the sub- 
array which results from feeding errors or mutual coupling between feedline and elements. The concept 
proposed in this paper was the design of 1x16 array from four 1x4 sub-arrays. The sub-arrays are 
manipulated in order to suppress the grating lobes by using parallel and series excitations. 
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2. PROPOSED SUB-ARRAYED PLANAR ARRAYS 

The array factor analysis of MxN (M=2, 4, 8, 16, 32) and (N=2, 4, 8, 16, 32) uniform broadside 
planar arrays which consists of an isotropic radiators are proposed in order to investigate the radiation 
characteristics of large-scale arrays. Figure 2 shows the conceptual diagram of describing this arrangement. 
Starting with 2x2 (4 element) planar array as a sub-array, a large planar array is formed up to 16x16 
(256 elements). For a planar uniform array of MxN elements in which all elements have same amplitude 
excitations, the normalized array factor is found as [3, 6-7]: 
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Figure 2. (a) 2x2, (b) 4x4, (c) 8x8, (d) 16x16 planar array arrangement 


The beam solid angle of planar array can be obtained as: 




( 2 ) 
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where 9 0 and 0 O specify the direction of the main beam, A 6 X is the half-power beam width (HPBW) of a 
broadside linear array where M number of elements and excitation of amplitude is same as the x-axis linear 
arrays forming the rectangular planar array while A Q y is the half-power beam width (HPBW) of a broadside 
linear array where N number of elements and excitation of amplitude is same as the y-axis linear arrays 
forming the planar array. 
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Directivity is defined as the ratio of maximum power density to its average value over a sphere. 
This shows how much the antenna is directive compared to isotropic antenna which has a directivity of equal 
to 1. To calculate the directivity of an array, the general expression shown in (4) is used [7]. 


Do 


= 4n 
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( 4 ) 


For large broadside planar arrays, (4) reduces to (5) 


D 0 = nD x D y cos9 0 


( 5 ) 


where, D x is the directivity of linear array placed in x-axis whereas D y is the directivity of the second linear 
array in y-axis. Moreover, the array solid beam angle kl A in (3) can also be used to find the approximate 
directivity of large broadside planar arrays as seen in (6). 

n 2 32400 

Un — — ( O ) 

^4[5r] n A [ de g 2] 


where kl A is mostly expressed in square degrees or square radians. The gain is calculated as: 

G = erD (7) 


So the gain is also defined as the ratio of intensity in a particular direction to the power gain of the reference 
antenna in its referenced direction. 


3. RADIATION PATTERN ANALYSIS OF SUB-ARRAYED PLANAR ARRAYS 

Based on the Array factor formula given in (1-7), a mat-lab code is created to plot the NXN (N=2, 4, 
8, 16, 32) 2D array factor plots starting with 2x2 (4 element) planar array as a sub-array and then a large 
planar array is formed up to 16x 16 (256 elements). The plots are as shown in Figure 3. 




(a) 




(b) 
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(c) 





(d) 

Figure 3. 2D, polar and linear plots of an antenna pattern with NXN planar arrays of isotropic 
radiators with a spacing of dx=dy=^/2, and equal amplitude and phase excitations 
a) 2x2 b) 4x4 c) 8x8 d) 16x16 


Figure 3 shows that starting from a unit array of 2x2, the pattern becomes more directive. The first 
unit array of 2x2 doesn’t show any side lobe but constructing 4x4 rectangular array, two side lobes suddenly 
appeared. Similarly, as we increase the modular array i.e. 2x2, 4x4, 8x8, 16x16, the number of sidelobes 
follow a specific pattern as 0, 2, 6, and 14 respectively. Therefore, the plotted array factor patterns clearly 
shows that the number of sidelobe for any NxN (2,4,8,16...) rectangular array can be easily predicted from 
the smallest 2x2 planar array. The combined for all planar array is shown in Figure 4. 



Figure 4. Combined array factor plots for 2x2, 4x4, 8x8, and 16x16 planar arrays 
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4. MODULARITY ANALYSIS OF SUB-ARRAYED PLANAR ARRAYS 

In this section, it is explained the results of research and at the same time is given the 
comprehensive discussion. Results can be presented in figures, graphs, tables and others that make the reader 
understand easily [3, 6]. The discussion can be made in several sub-chapters. 

4.1. Matlab analysis 

Based on the array factor characteristics presented in Figure 3 and using (1-7), the directivity, gain, 
3dB half-power beams widths, no of sidelobes and sidelobe levels are calculated to analyse the modularity of 
sub-array. The behaviour of modular planar arrays are closely observed to derive the impact of placing planar 
antenna arrays in modular form i.e. 2x2, 4x4, 8x8 etc. The directivity, 3dB half-power beam widths, number 
of sidelobes appearing and as well as the sidelobe level variations are presented in Table 1 and Figures 5 
and 6. 


Table 1. Comparison of directivity, gain, hpbw, number of sideslobes and sidelobe level variations 


No Array elements 

Directivity (dB) 

Gain (dB) 

3dB HPBW (degrees) 

No. of side lobes 

Side lobe level (dB) 

2x2 

7.06 

6.64 

60.33 

0 


4x4 

13.48 

12.67 

26.47 

2 

-11.3 

8x8 

19.71 

18.53 

12.87 

6 

-12.8 

16x16 

25.88 

24.33 

6.39 

14 

-13.15 

32x32 

32.11 

30.18 

3.19 

30 

-13.24 

64x64 

38.84 

6.51 

1.60 

62 

-13.26 



Figure 5. Variation of directivity, gain, no. of side lobes, and side lobe levels Vs no. of elements 



Figure 6. 3dB half-power beam width Vs no. of elements 
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Referring to Table 1, planar antenna arrays has shown similar repeating radiation characteristics 
when they are placed in a modular form. It can be observed some similarities in their directivities, half-power 
beam widths and sidelobes. To start with directivity, the directivity increases by 6 dB every time when the 
array size is increased by a factor of 4 such as 2x2, 4x4, 8x8, 16x16 and 32x32. For example, 2x2 planar 
array has a directivity of 7.06 dB while it has increased to 13.48 dB in the case of 4x4 array and so on. 
Similarly, the gain of the array antennas follows the directivity accordingly when an antenna with efficiency 
of 0.94 was used. Thus, it can be concluded that there is 6dB directivity and gain increase when the array 
enlarges along the given modular form. 

Moreover, the half power beam-width is exponentially decreasing with the increase of array size. 
The decrease in half-power beam width is indirectly revealing the increase of that of directivity as the array 
becomes larger. The 3dB HPBW variation shows a halving trend where it drops half every time when array 
size is increased modularly. For instance, 3dB HPBW’s for 2x2, 4x4, 8x8 planar arrays is 60.33, 26.47, and 
12.87 degrees respectively. 

From Figure 5, it is clearly seen that the directivity and gain are linearly increasing with the increase 
of number of elements in the array. Similarly, sidelobe level variations also increase linearly in the beginning 
before it becomes constant. But the number of sidelobes and 3dB half power beam widths follow an 
exponential form where the directivity and HPBW are inversely related as illustrated in Figure 6. 

4.2. CST analysis 

Previously, Matlab software was used to analyse the radiation characteristics of a rectangular 
isotropic planar arrays to observe the effect of developing large arrays by repeating a small rectangular sub- 
array many times placed in a modular form i.e. 2x2, 4x4, 8x8, 16x16 etc. It was found that those isotropic 
planar arrays has shown similar repeating radiation characteristics in their gain, directivities, half-power 
beam widths and sidelobes. Moreover, a mathematic model was derived to predict the radiation performance 
of large-scale arrays from the small sub-array. 

Therefore, in order to validate those findings, CST software package is used to simulate planar of 
patch antennas for validation. In this part, the practical antenna effects such as mutual coupling between 
elements will be considered. Therefore, the results from the Matlab and CST will be compared in order to 
validate the derived model to predict the radiation performance of large arrays such as array factors, 
directivities, gains, half power beam-widths, sidelobe ratios and so on. 

A single rectangular array was constructed and showed in Figure 7. The antenna is designed using 
Rogers RT5880 (lossy) dielectric substrate having with a dielectric constant is 2.2 and a thickness of 1 mm. 
The feed line is a 50 D microstrip line attached to the patch. The substrate is fully grounded. The operating 
frequency is set to 6 GHz. The behaviour of planar arrays are observed in 2x2, 4x4 and 8x8 on some specific 
array characteristics such as the directivity, gain, 3dB half-power beam widths, number of sidelobes 
appearing and as well as the sidelobe level variations and Figure 8 represents the structure of planar arrays. 
The results are compared in tabulated form as shown in Table 2. hus, this analysis can be utilized to predict 
the radiation characteristics of large-scale arrays. 



0 



(a) (b) 

Figure 7. Single planar array results: (a) Structure, (b) Pattern 
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S-Parameters [Magnitude in dB] 
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(c) 

Figure 7. Single planar array results: (c) Sparameter 
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Figure 8. Structure of planar arrays, (a) 2x2, (b) 4x4, (c) 8x8 


Referring to Table 2, after comparing theoritical calculation results and simulation results it can be 
found some connections between them. In terms of directivity the simulated results are minimum 4dB higher 
than the calculated results. Moreover in the simulation results for 2x2, 4x4 and 8x8 almost 6dB increase of 
directivity can be observed. In terms of gain the simulated results are also higher than the calculated results. 
For 3dB HPBW the degrees drop in calculated results are much drastric than the simulated results. However, 
number of sidelobes are same in calculated resuls and simulated results. Table 3 illustrate 2D radiation 
pattern of theoretical calculation and simulation results. The similatiries between theoretical and simulation 
are clearly visible as the main lobes of 2x2, 4x4 and 8x8 are more likely similar and shoed the similar 
number of sidelobes. 


Table 2. Comparison of directivity, realized gain, 3dB angular width, 
number of sideslobes and Sidelobe Fevel Variations (CST) 


Array 

element 

s 

Directivity (dB) 

Realized Gain (dB) 

HPBW (degrees) 

No. of sidelobes 

sidelobe level (dB) 

Calculate 

d 

(Matlab) 

Simulate 

d 

(CST) 

Calculate 

d 

(Matlab) 

Simulate 

d 

(CST) 

Calculate 

d 

(Matlab) 

Simulate 

d 

(CST) 

Calculate 

d 

(Matlab) 

Simulate 

d 

(CST) 

Calculate 

d 

(Matlab) 

Simulate 

d 

(CST) 

2x2 

7.06 

12.6 

6.64 

8.94 

60.33 

29.8 

0 

0 


-13.7 

4x4 

13.48 

17.7 

12.67 

13.7 

26.47 

23.4 

2 

2 

-11.3 

-14.6 

8x8 

19.71 

23.1 

18.53 

19 

12.87 

13.6 

6 

6 

-12.8 

-13.2 
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Table 3. Comparison in term of 2D radiation pattern between theoritical calculation and simulation 

Simulation- CST 
Single patch antenna 


Specification 


Theoritical- Matlab 


Simulation-CST 



5. CONCLUSION 

Large-scale antenna arrays have been concentrated to assess the possibility of a basic and productive 
method of utilizing sub-arrays for the improvement of expansive uniform arrays. Large arrays can be 
fashioned if small sub-arrays are repeated throughout the remaining large array. The antenna sub-array is 
examined based on the pattern, the directivity and the gain, 3dB half-power beam width, number of side- 
lobes and also the side-lobe level. There is 6dB increase in gain and directivity when the array is enlarged 
along with the given modular form. The modularity is maintained in the CST simulations of the sub-arrays 
and validate the matlab analysis. By increasing the size of the array the half power beam-width decreases 
exponentially. The variations caused by the level of the side-lobes also escalates linearly at first before 
becoming constant. But the quantity of side-lobes and 3dB half power beam widths follow an exponential 
method. So it can be concluded that, this analysis can be used to calculate the radiation characteristics of 
large-scale arrays by forming a small sub-array of 2x2 components as a perspective. 
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